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1 Introduction. 

Statistical approaches to the equations of fluid dynamics are well known 
[10,11]. The present author, following a different stochastic approach to the 
previously quoted classical works, has given random implicit representations for 
the Navier-Stokes equations for an incompressible fluid on a smooth compact 
manifold without boundary. The approach chosen for the obtention of these rep- 
resentations stemmed from Stochastic Differential Geometry, i.e. the invariant 
theory of diffusion processes on smooth manifolds, and in particular, in Euclidean 
space [1,2,8,9,14,23-27]; this approach has lead the author to give as well random 
representations for the kinematic dynamo problem, i.e. the equations of trans- 
port of a passive magnetic field transported by the fluid, in the same geometric 
situation [5]. While stemming from the fussion of gauge theory and stochastic 
processes on manifolds, the present theory is founded on the relation between 
a class of linear connections with torsion [13,23,27], and the random diffeomor- 
phisms determined by them, extending thus the classical diffeomorphims associ- 
ated to perfect fluids, presented in the pioneering work of V.I. Arnold [15,16] and 
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further ellaborated by Ebin and Marsden [17]. In a recent interesting program, 
Gliklikh [18-20] has extended the original proposal of Arnold-Ebin-Marsden, to 
give a formulation of the Navier-Stokes equations on the n-torus as a stochastic 
perturbation of the group of volume preserving diffeomorphims associated to the 
flow of perfect flows. In the present gauge-theoretical approach, the representa- 
tions obtained by the present author were obtained in two instances [5,6,12,13]: 
Firstly, using the so-called derived (or jacobian) random flow, constructed as the 
derivative flow of the lagrangian random representations for the fluid particles, 
in case the manifold is isometrically immersed in Euclidean space (the case of 
gradient diffusion processes). Secondly, for the more general case of arbitrary 
compact manifolds, the representations were derived in terms of the Hessian and 
Ricci random flows related to the former processes. In this note we wish to extend 
these constructions of the random representations for the vorticity Navier-Stokes 
and kinematic dynamo equations, to the case of a smooth manifold with smooth 
boundary, both for which no general representations are known; in this case, the 
boundary conditions for the vorticity in practice in Fluid-dynamics, will turn 
to be the well known reflection at the boundary conditions which we shall rein- 
troduce below. Yet, our presentation will follow instead the original approach 
to Brownian motion on smooth manifolds which lead to Stochastic Difi^erential 
Geometry: The stochastic extension of the Cartan classical soldering method, 
as originally presented by Malliavin [8] and Elworthy [2], which we have pre- 
sented as the mathematical approach to Statistical Thermodynamics, Quantum 
Mechanics and Gravitation [4,13,23-27,29], while the classical method appeared 
to be adequate for the formulation of the dynamics of classical relativistic spin- 
ning particles subjected to gravitational fields [22,28] without using no lagrangian 
nor hamiltonean structures. In this setting, the torsion is no longer provided by 
the full irreducible expression of the torsion tensor, instead, the trace-torsion is 
the basic 'physical' field to be taken in account in a theory of generalized Brow- 
nian motions, since its conjugate vector-field turns to be describe the drift of the 
process. Yet, while in the setting of classical differential geometry in which the 
transport of fields is done along smooth curves so that the Cartan calculus of 
differential forms is the basic tool, in this continuous non-differentiable setting 
this calculus has to be supplemented by the rules of stochastic analysis. These 
together yield a most formidable theoretical and computational instrument, and 
the Ito formula for differential forms, its highest and most basic expression of this 
extended calculus. Thus, it will turn to be of no surprise, that it is precisely this 
formula which will give us the key to the derivation of the analytical expression 
for the evolution Navier-Stokes equations for the vorticity on a manifold with 
smooth boundary. This approach will turn out to extend not only the analytical 
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expressions obtained in the case of no boundary case obtained by the author, but 
the original approach in Computational Fluid-dynamics known as the random 
vortex method, which deals with viscous fluids on the Euclidean plane [40]. 

2 Riemann-Cartan-Weyl Connections and Their Lapla- 
cians 

We shall assume an n-dimensional smooth connected oriented compact man- 
ifold, M. Wc shall further assume as given a Riemann-Cartan-Weyl connec- 
tion on M, V, i.e. a linear connection on M (or still, a covariant derivative on 
TM) which is ^-compatible, and such that its Christoffel coefficients defined by 
V_^gfT = V^^^ are of the form 

r^^(ar) = I + {6^ Q^{x) - gp^ix) Q'^(x)} ,x G M (1) 

with torsion tensor 

T|,(x) = i/2(r^^(x)-r^^(x)), 

which has an irreducible decomposition which is zero for all its components with 
exception of the trace-torsion 1-form given by 

Q{x) = Q0ix)dx^ := T^p{x)dx^. 

The first term in (1) stands for the Christoffel coefficients of the unique torsion- 
less linear connection determined a Riemannian metric g, i.e. the Levi-Civita 
connection, which we shall denote in the following as V^. 

Associated to this connection we have a Laplacian operator defined by 

i/o(5,Q) :=1/2V2 (2) 

(the subscript denotes its actions on functions, which wc here consider to be 
smooth, as well as all other fields that appear subsequently) which only depends 
on g, its derivatives, and Q, since one can compute explicitly to obtain 

Ho{g,Q) = l/2Ag + LQ,. (3) 

In this expression 

Ag = div grad = (V^)^ 
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is the Laplace-Beltrami operator associated to V^. Furthermore, is the vec- 
torfield conjugate to Q by g, i.e.Q^{f) = g{Q, df), for any smooth function /, and 
Lx denotes the Lie-derivative with respect to the vectorfield X on M. Recah 
that on functions / we have the identity Lggf = Q^{f), so that by recalhng that 
the metric is non-degenerate, we conclude that expression (3) corresponds to that 
of the most general laplacian with zero-potential term, acting on functions. We 
note that the second term in expression (1) describes the drift, and appears that 
in an invariant setting, it is related to the trace-torsion 1-form. 

These operators can be rewritten in a form suitable to carry out stochastic 
analysis on M and on Po{n)^ the bundle of orthogonal frames of the tangent space 
TM over M [1,30], i.e. 

d 

-fb(n) = {r = {x,e{x)),x G M,e{x) = e^(x)^^|^ a basis of T^M} 



so that 



jeaCf = Sab- 



Consider the canonical horizontal vectorfields on -Po(n)> ^^'^ defined by 
the connections V and respectively, with components L'^ ,L^'' ,a = 1, . . . ,n 
respectively, given by Q 

L^Fir) = e:{x) — ir) - r^^(x)e2(x)ef (x) — (r), a = 1, . . . ,n (4) 

and 

LrF(r)=e^(x)^(r)-|jJ(x)e2(:r)ef(:r)^(r),a = l,...,n. (5) 

Then, the horizontal lift of the operator Ho{g, Q) acting on functions defined on 
-fb(n)) which we denote as HQ{g,Q) is (see Theorem 1.2, page 238 in [1]) 

Ho{9,Q) = Il^{L^), (6) 

so that for any "basic" function, /(r) = / o 7r(r) = f{x),r = {x,e), where 
TT : Po(n) ^ ^(^) '■= X foi^ any r = (x,e) € Po{n) is the bundle projection, we 
have the identity 

[Ho{9.Q)f]{r) = [H^{g,Q)f][x) (7) 



^Einstein repeated sign sum convention is assumed in the following. 
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and the horizontal hft of the Laplace-Beltrami operator, consequently verifies 

A, = 2Ho{9,0)=L^'{Lj'). (8) 

Note that we also have 

Ho{g,Q) = l{L^'{L^')) + LQ„ (9) 

where Q» is the horizontal lift of QS, i.e. QS(/) := (Q^(/)) = vr*(Q''(/)). 

We can now extend these operators to act on smooth differential forms de- 
fined on M and its horizontal lifts to act on smooth differential forms on -Po(n)- 
Indeed, consider 

Hp{g, Q) = -Ap + Lqg | ap(t*m) , (10) 

where 

Ap = (d - 6)'^\{KP{T*M)) 

is the Hodge laplacian acting on differential ;?-forms defined on M, so that d and 
5 are the exterior differential and co differential operators, respectively. We recall 
that 5 is the adjoint operator to d with respect to the pairing 

where volg{x) = det {g)dx^ A ... A dx'"" is the Riemannian volume density, (/> = 
^'Pai...ap{x)dx'^^ A ... A dx'^p and u = ^i^aiA...Aap{x)dx°'^ A ... A dx^p are p- 
forms and = gl^f^^ ■ ■ ■ g^^"'^(t>i3i...i3p- Since (f = 0, then also (5^ = and 

consequently 

Ap = -{d6 + 6d)\/^p(T*M) (11) 

and then 

Hp{g,Q) = ~{6d + d6) + LQg\f^p(^J^,M)■ 
As well known, the Lie-derivative of a p-difi'erential form with respect to a vec- 
torfield is well defined and independant of the metric g: 

Lx<^ = {ixd + dix)u) 
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where ix is the interior derivative with respect to the vector field X and lo is an 
arbitrary p-form on M , and in case of p = 0, since ixf = for an arbitrary scalar 
field / {ix reduces degree) we get the original laplacian in (3). While the second 
term is independant of the metric, the Hodge laplacian depends on the metric 
and the curvature obtained by taking appropiatc expressions on its derivatives. 
For the forthcoming discussion it is essential we display the dependance of the 
Hodge laplacian on the curvature, to write the whole laplacian of our interest 
given by (10) in the form 

1 1 ^ 



- E (-l)'+''^^a::a,<^7/3ai...d,...d....a, + (i^Q«</')ai...a,. (12) 

l<^<i'<p 

Here, as usual, a denotes ommision of the index a. The second and third terms in 
the r.h.s. of (12), constitute the well known Weitzenbock term, which of course, 
vanishes completely in the case p = 0. Note that for 2-forms, 

— 2" ^^CKl CK2 dx /\ dx ^ 

the Weitenzbock term is 

- \{-^r E - Y.^-ir^'^i^^i^^^.p (is) 

i/=i ^l<v 

where in the first term we have a coupling of the differential form to the Ricci 
curvature tensor, Ric = {R^) = (R^-y'Ja), and R = (-R"'^') is the (2,2) tensor 
associated to the Riemannian curvature tensor [31]. 

Now, given a j>form defined on M, we can define a p-form on Po(n)' its 
horizontal lift, and thus further establish an isomorphism between Af{T*M) and 
hor {AP{T* Po(n)))i the horizontal p-forms defined on -Po(n)) i-6- the subspace 
of AP{T* Pqi^^)) whose projection by the bundle mapping tt has zero kernel, i.e. 
7r*(hor {APT*Po(n))) = A^(r*M). Indeed, consider the co-frame bundle P^^^^ = 
{{x,9{x) = (6'"(x)) = (6'2(x)o?x° G T*M),9''{x){ef,){x) = (5^, for any x G M}, so 
that eie^pg^'l^ = 6ab- Then, given a p-form n on M, 

n = n{x) = -^na,...ap{x)dx'^' A ... A dx'^", (14) 
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we can define functions on Po(n) given by 

^a,...a,ir) := na,...a,{x)e2l{x) . . . e-;(x),r = (x,e) (15) 

and thus obtain a (horizontal) p-form on Po(n)) i-^- an element of hor(A^'(r*Po(„))) 
defined by 

n = n{r) = lj7„,...„^(r)ri A ... A (16) 

This definition can be inverted: if we start with Cl G APr*(PQ(„)), we can obtain 
the p-form fl defined on M by 

n = n{x) = ^na,...ap{x)dx'^' A ... A dx""", (17) 

where 

^a,...a,ix) = ^a,...a^r)ei\[x)...eX{x),T = {x,e),e = e-\ 

Naturally, defined by (15) is 0(n)-equivariant, i.e. 

n„,..„,(r) = (l,,...,^{TAr)Bl\ . . . (18) 

for any r G Po{n)i ^ & Po{n)^ B the inverse of A and T^r = Ar. 

Suppose given a p-form Cj defined on Po{n)'i then, we can consider a time- 
dependant p-form J7 defined on [0, co) x Po{n)i by the expression 

^ai...a^{T,r) = Er[u:ai,„a^{r{T,r,w))], (19) 

where E denotes expectation value, w : [0, oo) — > i?" is a continuous mapping such 
that w{0) = is the canonical realization of a Wiener process, and r = r(T,r,w) 
satisfies the Stratonovich stochastic differential equation 

dr{T) = L^{r{T)) o dw{T), r(0) = r. (20) 

We note that if r(0) G -Po(n)) then r{T,r,w) G Po(n))^''" ^ 0- Thus defined, 
the scalar components ^ai...ap of CI{t, r) = ^^ai-apO""^ A . . . A9"-p e Ap([0, oo) x 
P*Po(n))i < ai < . . . < ttp < p)aie the unique solution of the Cauchy problem 
(see page 286 in [1]) 

^ = Ho{g,Q)V (21) 
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with 



V{0, -) = VI < ai < . . . < Op < p. (22) 

We wish to solve, more generally, the Cauchy problem for p-forms on [0, oo) x 

Po{n)- 

— {T,r) = Hp{g,Q)n{T,r), (23) 

with given 

niO,-)=cjeAPT*Po^n), (24) 

where we notice that in distinction with the problem (21&;22), a coupling of the 
Weitzenbock term to the components of O is to be accounted. This problem is 
simply solved through the use of the Feynman-Kac formula. For this we consider 
the canonical realization of the generalized horizontal Brownian motion, given by 
the set 



T^(Po(n)):=Co([0,oo)^Po(n)) 
of all continuous mappings 

w : [0,oo) Po(n)> 

where B(VF(P(5(„))) is the cr-field on M^(Po(n))) generated by the Borel cylinders 

and Br(W{Po(n))) be generated by the Borel cylinders up to time r. Let Pr be 
the probability space law on W{Po(n)) of r ^ r{T,r,w). Consider for each w G 

W(Po(n)), the time-dependant (j),p)-tensoic M{t,w) = {Mal'.'.'.apiT,w)) defined 
on [0, oo) X Po(n) as the solution of the system of equations 



dMt.:l, ... 1 



^{w{r)) = -^(-l)^^Ric^^^^(r))M„^^^if ''+--^''+^(.i;(r)) 



dr - T 



- J2 (-l)"+''^'1'(^(r))M„'^-:i+^^'''+^-''''''''+^-'^(u;(T)), (25) 
i</i<i/<p 

with 

(0) = S',\...6^,l,yi<ai<...<ap<n, (26) 
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and R := (^"^d) = (-^°^//e^e^6'26'^), the horizontal Hft of the Riemannian curva- 
ture (2, 2) tensor, which we shall write as R = R®e®e®6®9, so that we also 
have that R = R® 6 9 ® e ® e, axid Ric^ = (Ric^ \) = {RID = {R^ef^el) which 
we shall write as Ric'' = Ric ® e® 9,is the horizontal lift of the Ricci tensor, to 
Po{n)- For later use, we introduce the tensor 

so that 

R^®e = R. 

It is clear that system (25, 26) has a unique solution. Prom the uniqueness 
of the solutions, it follows that 

= (27) 

for every permutation a. One can prove still that this solution is 0(n)-equi variant, 
since forall A G 0{n),w G W{Po{n))i 

MIWX ^^^) = ■ ■ ■ bXaI\ . . . AX , (28) 

where B = A~^. Let now 

io = uir) = A . . . A e''- G (r*Po(„)) 

be the horizontal lift of a given lo = lij{x) G A^(T*M), and consider the differential 
form U, defined on [0, oo) x Po(n) with components Uai...ap{''',r) given by 

U{t, r) := Er[M{T, w{t)) ® oj{w{t))] (29) 

which componentwise is 

Ua,...a,ir,r) = S,Ki;;.i^^(r,u;))i;,,...,^(u;(r))]. (30) 

This expression is alternate in (ai, . . . ,ap) and obviously, it is 0(n)-equivariant. 
Consequently, there is a p-form, U, defined on [0, oo) x M, whose horizontal lift 
is U. If we consider the time dependant p-form given by 

M®{v:= ^Mt^:ilOb,...bX' A . . . A r^ 
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and apply to it the Ito formula, we get 



(M (g) tj)(r, w{t)) - (M (g) c;;)(0, r) = 
a martingale + / M [— — \- Hp{g,Q)Lb)\{s,w{s))ds. (31) 

Therefore, 

U = U{r, r) = ^Ua^.-a^O"' A ... A ^-(r, r) 

with components given by (30) is the unique smooth solution of the Cauchy 
problem 

^ = Hp{g,Q)U, (32) 

with initial condition 

C7(0,r) =u;(r). (33) 

Theorem 1. The p-form U G A*'([0, oo) x T*M), such that its horizontal 
lift coincides with U, i.e. U = (U{Dtt)) = U fS>'^e, where D denotes the derivative 
and (g^e denotes the p-th tensor product with e, is the unique smooth solution 
of the Cauchy problem 

— = Hp{g,Q)u,uiO,-)=u;. (34) 

This theorem extends a construction originally due to P.Malliavin, for the 
Riemannian connection driftless case [1,7,30]. The process M constructed above, 
to give account of the Riemannian curvature terms, is called a multiplicative 
operator functional. 

To conclude with these definitions, we note that from (6, 10&:12) follows an 
expression -which we shall use later- for the laplacian on 2-forms (j){r) = ^a?,^" 
defined on Po(n) '■ 

Hp{g, Q)(t>ab = ^ (L^)2</,„, + i (RicV)a6 - (i?'^)a&- (35) 
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3 The heat equation for manifolds with boundary 



Now, we shall assume further conditions on M, that of having a smooth 
boundary dM. Near the boundary, we can choose a coordinate open neighbor- 
hood U of M, and smooth coordinate functions (x^, . . . such that > 0, 
for all X G [/; furthermore x G ?7 fl dM if and only if x"^ = 0, and we shall further 
assume that we have a smooth Riemannian metric g = {gap) on T*M such that 
Qanix) = 0, for a = 1, . . . ,n - 1. 

Let be a p-form on M, i.e. ^ G A*'(r*M). Denote by 0tan = <i>\dM-, the 
tangent component of (j). and the normal component by (/>norm = 'p — 4'tim- We can 
still represent the normal component using duality. Indeed, if * : AP{T*M) — > 
A"'~^'(r*M) is the star operator defined by the metric 5, then 

(_l)p(n-p)+p-l[^(^^ A r?) A ry]|aM = <Anorm, (36) 

where r]{x) is a 1-form such that 'q\dM = gap{x)Ti^{x)dx°',x G 9M, where 

n(x) :=n"(x)— = - — ,a = l,...,n,x G 5M, (37) 

is the inward pointing normal unit vector to dM at x. 

Then, we shall say that a p-form is said to satisfy the absolute boundary 
condition if 

'?!'norm = 0, (d^)norm = 0. (38) 

Note that in the case of a 2-form on M, $7 = ^U^fidx'^ A dx^, the absolute 
boundary conditions for O, are 

Qan = 0,a = l,...,n-1, (39) 

and 



ax" 



0,a,/3 = l, ...,n-l, (40) 



respectively. 

We wish to solve the boundary-initial-value problem for 

^ = l^a^aAr,r)9''' A Oa, G A2([0, Oc) X r*Po(n)) 

given by 

— (r,r) = i?p(5,Q)J7, n{0,r) = u;{r),r e Pq^u), (41) 
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where a) is the horizontal hft to -Po(n) of ^ given p-form, cv defined on M, with 
the absolute boundary conditions. In fact, conditions (39&;40) can be lifted to 
[0, oo) X Po(n) to give: 

W^^aia.(r,r)|aPo(„) =0,e = e-' = {9-Jx)dx-), (42) 



C,(^)C.(^)^f^a,a.(r,r)|aPo(„) = 0,a = 1, . . . ,n - 1. (43) 

where dPo(n) = {r = (x, e) : x G dM}. We shall solve problem (41, 42&43) using 
the generalized horizontal Brownian motion on Po(n) a multiplicative hori- 
zontal operator functional. We start by reformulating the boundary conditions. 
Let P = (P^) defined by P" = 1, and Pp = 0, for a,(3 ^ n; consider further 
Q := I — P. Thus, the absolute boundary conditions (42&;43) can be rewritten 
as the single equation 

P ^ e (g) e (g) n{T,r) + Q e (g) e (g) —{T,r) = 0,r = {x,e) G dPo{n),G = e~^(44) 

Furthermore, we can rewrite the boundary conditions (44), instead of as a con- 
dition on A^([0, oo) X r*Po(„)), as a condition for 

Q.\T,r) = Qt(r) = ^Q^„(r, r)dx° A ^(x) € L{TM,R''*),W > 0, 

i.e. for each r > we have a linear mapping from TM K^* such that for each 
r = {x,e) G Po(n), ^rir) maps T^M into hor (r*Po(„)) ~ R"^' , related to j^(r,r) 

by 

e<gn^ = n 



and satisfying the boundary conditions 



P (8) n',{r) + Q®e®[—^{r)- ^(r)r^^(x)e^(ar) 



m 



+ e{x)^rn{x)^9{x)®ni]{r) =0,Vr = (x,e) G dPo(n), (45) 

where r„(x) = (r^^(a;)), and (9 ® ^^)a/3 = (^^/j^^aa) multiplication in the 

representation space P"* of this linear mapping. Indeed, note that if Ct{T,r) is 
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0(n)-equivariant, then 0''(r, r) = J7(t, x) (8>e(x),r = (x, e) E -Po(n) Q> and thus 

^r^^(x)e^(x) = e®r„,(x) 0(g) Ot- (46) 

Therefore , in account of identity (46), by further multiplying (45) by 9, we 
retrieve the boundary conditions (44) for equivariant 2-forms in A^([0,cxd) x 
r*Po(-„)), or multiply by 6 on the right, to obtain a 2-form in A^([0, oo) x T*M) 
as explained in the previous footnote. 

From a localization argument, we can assume for simplicity that M is the 
upper half-space of i?", {x : x = (x^, . . . , x") G R"-} and dM = {x G M : x*^ = 0}. 
The metric g on M is assumed further to be bounded and smooth together with 
all its derivatives, i.e. G C^(M)., 1 < a < 13 < n. 

We start by considering the following Stratonovich stochastic differential 
equation for the diffusion process r(r) = (X(r),e(T)) on Po(n) ~ -?^+ -R" i 
following [1]: 

dr°(r) = L^(r(r)) o di?'=(r) + 5°#(r), r(0) = r, a = 1, . . . , n, (47) 

where as before, V denotes a Riemann-Cartan-Weyl connection with coefficients 
of the form (1), L^, k = 1, . . . ,n denote the canonical horizontal vectorfields de- 
fined by this connection on Po(„) in (4), and B denotes a n-dimensional Brownian 
motion. We note that it can be rewritten, for r = {X, e), as the system 

dX"(r) = e^(r) o dB^{T) + (5°#(r), X(0) = x, a = 1, . . . , n (48) 



det{T)= - T^p^{X{T))el{r)odXP{T) = 

- T%{X{T))el{T)el{T)odB^{T) 

- r°^(X(r))e2(r)#(r), e"(0) = e(x) G T^M, a. A; = 1, ... , n(49) 

Here, i;A(r) is a continuous non-decreasing process which increases only when 
X{t) G dM, i.e. a local time on 5M, and which causes the reflection of the 
process at the boundary. 

This is a generalized (i.e. with non-zero drift described by (Q^), the hor- 
izontal lift of Q^) horizontal Brownian motion on the bundle of orthonormal 
frames Po{n) with reflecting boundary. Indeed, it can be easily seen that if 

^This identity, which can be rewritten as f2(r, x) = (r, r)^6(x), or still f2(r, r) = f2(r, x) ® 
e ® e, will play a crucial role for obtaining the expression for the vorticity as defined on M 
instead of Po{n))', we shall work out the boundary-initial- value problem for Navier-Stokes on 
[0, oo) X L{TM,T*Po(„)), to later project on A^([0, oo) x T*M) by this procedure. 
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(X(0),r(0)) € Po(n), then (X(T),r(T)) G Po(n)- We readily note that if 9M = 
and if we further set = in eq. (47), we obtain the generaUzed horizontal 
Brownian motion on Po(n)) for a manifold without boundary. Having assumed 
that gafs G C^{M), if we assume further that Qa e C^{M) (1 < a < ra) and 
(f) G C^{M), then eqt. (47) has a unique strong solution (c.f. Theorem 7.2, 
chapter IV, [1]). 
Set now 

A--{r) = 5""(x), (r) = -e'l{x)T-p{x)g--{x),r = (x,e). (50) 

The following result is crucial to our following constructions. 

Theorem 2. Let r(r) = (X(r),e(T)), the generalized horizontal Brownian 
motion with reflecting boundary condition described by eqt. (47). For any 
smooth function F{t, r) defined on [0, oo) x Po(n) we have 

dF(r,r(r)) = [{Ho{g,Q)F){T,r{T)) + ^{T,r{T))]dT 

+ {L^F){t, r{T))dB\T) + (X„F)(r, r(r))#(r), (51) 

where X„ is the canonical horizontal lift of Xn = i.e. 

dF A"'"(r) dF 

(^nF)(.,.) = ^(.,0 + ^g^(r.r). (52) 

Proof. It follows from the Ito formula that 

dF(r,r(r)) = e^(r)^(r, r(r)) o di?"(r) 

- ^(r,r(r))r^,(X(r))eJ(r)ef(r)od5»(r) 

dF dF 
+ -^(r,r(r))dr + — (r,r(r))#(T) 

- ^^{r,r{r))KAr{r))e:dct>{T) 

= (LjF)(T,r(T))odi?'^(r) + — (r,r(r))dT + X„(F)(r,r(r))#(r) 

so that 

d{L^F){T,r{T)) = LfiL^F){r,r{T))odB\T) 

+ ^^^(^' ^(^))^^ + Xn{L^F){r, r{T))d4>{T) 
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and then 

(LjF)(T,r(T) odB(T) = +^diL^F)iT,riT)).dB'^iT) + (LjF))r, r(r))di?"(r) 

= l{L^)'F{T,r{T))dT + {L^F){T,r{T)) . dB^r) 

= ^Ho{9, Q)F{t, r{r))dT + {L^F){t, r{r))dB-{r), 

the last identity follows from (6), and • denotes the Ito contraction on stochastic 
differentials [1] , with which we can conclude with the proof. 

Prom Theorem 2 follows that the process r(T) G -Po(n) is determined by the 
differential generator Hq^q, Q) with the boundary condition X^F = on dPQ^^y 
Furthermore, 

dX"(r) . dX^'ir) = ^"'*(r(r))dr, (53) 

dX"(r) • delir) = A'^''"(r(T))(ir, a, = 1, . . . , n. (54) 

To obtain a solution of the boundary-initial-value problem (41&45) for the 
case of smooth boundary, we need to construct a multiplicative operator func- 
tional, which extends the one presented for the boundaryless case presented 
in Theorem 1. Consider a canonical realization of the generalized horizontal 
Brownian motion on Po{n) with reflecting boundary condition. Wc thus con- 
sider W{Po{n)) = ^"([0,00) Po{n))j provided with the probability measure 
Pr of the solution of (46) and, r{T,w) = w{T),yw G W{Po(n))- Given a Borel 
probability measure /x on Po(n)) we define Pfx{B) = Jp^^ ^ Pr{B)iJ.{dr),yB G 

B{W{Po(n)))- Let = n^BiWiPo^n)))^" and T,, = {A £ F: for any /i, there 
exists 5^ G Br{W{Po(n))), such that P^{AAB) = 0}. In the following we fix 
fi and consider the probability space {W{Po(n))i^iPiJ,)- Writing for simplicity 
r(r) =r{T,w) = {X{T,w),e{T,w)) = (X(r),e(T)), we set 

0(r) = lim^^oY^ £ I^o,e){X^{sW{x)ds, (55) 

and 

B\t) = r e;(s) o [dX^is) - Sycl^is)], (56) 
Jo 

where 6 = {9^{x)dx°') is the dual (co-tetrad) field to e. Then, {B^{t) : i = 
1, . . . , n} is an n-dimensional (^T-)-Brownian motion and ({r'(r) = r(T, w)}, 4>{t)) 
satisfies the s.d.e. equation (47). 



15 



Lemma 1. {Pr} is invariant under the action of 0(n), i.e. if we define 

{Taw){t) = TA{w{r)),\/w G w;(Po(n)),T > (57) 
and if TA{Pr) is the image of Pr under the mapping w —>■ Taw, then, 

TaiPr) = PtMv)- (58) 

Proof: Let r(r) be a sohition of eq. (47) with B(t) and 0(t) such that 
r(0) = r. Then for A G 0{n), f{T) := TA{r{T)) is a solution of eq. (47) with 
B{t) = A~^B{t) and 4>{t) = 0(r), such that f(0) = Ar. If we further observe 
that B{t) is another n-dimensional Brownian motion, then TA{Pr) = PxAir) holds 
by the uniquenesss of solutions. 

As a consequence of Lemma 1, we note that X{t) = 7r(r(r)) defines a diffu- 
sion process on M whose differential generator is Ho{g,Q), and with boundary 
condition given by the vanishing of the normal derivative of / at x G dM: 

df df 

-± := „"(a;)-^ = 0, X G dM. (59) 

This diffusion is the so-called generalized Brownian motion on M with reflecting 
boundary condition; it represents the random process on actual configuration 
space, accounting with the given boundary conditions. 

In the following we shall fix a probability measure /x and consider the prob- 
ability space {W{Po(n)),^,PfM)- Then r(r) = r{T,w) = {X{T,w),e{T,w)) = 
(X(r),e(r)) is the solution of eq. (47) with B{t) and 4>{t) given by (55&56). 

We shall consider the i?" -valued process K(t) = {K^(t,w)) defined 
-in extending a construction due to H. Airault (see [3])- as follows: Set K^{t) = 
K{t) ® P and K'^{t) = K{t) ® Q; hence K{t) = K^{t) + K^{t), where 

(z)for any r > such that X{t) G int(M), 

dK^{T): = dK{T)®P = K{T)®[d{T)(^de{T) 

+ {^Ric{X{T))-R{X{T))}dT]®P, (60) 

(M)for any r > , dK'^{T) := dK^{T)Q = K{t) ® [^(t) de(r) 

+ {^Ric(X(T)) - i?(X(r))}dr]/int(M) {X{t)) ® Q, (61) 

and with probability 1, t —>■ K^{t) is right-continuous with left-hand limits; 
furthermore 

{Hi) K^{t) = 0, ff X{t) G dM, (62) 
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and the initial value is 

(iv) i^HO) =/int(M)(^(0))e(0) 55P,K2(0) = e(O)0Q. (63) 

Note that in terms of a component representation for K, we have that the above 
expressions appear to be of the form 

idK)i^P^ = K^[i9^de)}P^ + {R} + RjJPs^dr 

= K^[i9 de)lP^ + {e^Rtel + e'^R'^leDP^dr]. 

and similarly for K^. |^ 

An equivalent formulation is the following. An i?" ® i?" -valued process 
adapted to {J-t) is a solution of the above s.d.e. with the given initial conditions 
iff 

K^{t) = K{T)®P=h^^M^)+ f K{s)®[e{s)®de{s) 

Jo 

+ {^mciX{s))-R{X{s))}ds])^P 
+ /|,>^| r K{s) [9{s) (g> de{s) + { jRic(X(s)) - R{X{s))}ds] O P, (64) 

~ Jt(T) 2 

K^{t) = K{t) ®Q = e(0) 0Q+ r K{s) ® [9{s) de{s) 

Jo 

+ {^Ric(X(5)) - RiXis))}ds]I,^,^M)iXis)) Q, (65) 

where 

a = inf{s : X(s) G 9M}, (66) 

and equal to if this set is empty, is the first-hitting time of {X{t) : r > 0} to 
dM, and 

t{T) = sup{s : s < r, X{s) G dM}, (67) 
and equal to if this set is empty, is the last-exit time from dM. 



It may seem that instead of taking R we should take the tensor with components {R^^] 



(R^s^d^j yet on muhiplying with P and Q it turns to be indistinct with our previous choice. 
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Let E be the set of all K^^R"' -valued processes ^(r) defined on (VF(Po(n))) 
P) adapted to ( J-t) such that r — > ^(r) is right continuous with left-hand 
limits almost surely and satisfies 

sup,g[o,T]^M[||C(r)||'] < oo, VT > 0. (68) 
Define a mapping $ : S — > S 

$He)(T) : = mir)P = I{r<aMO) + £ Kis) [d{s) ® de{s) 
+ {^Ric{X{s)) - R{X{s)))ds]) ® P 

+ /|^>^. K{s) ® [e{s) ® de{s) + (^Ric(X(s)) - R{X{s)))ds\ (g) P, (69) 

~ Jt(T) ^ 

$2(^)(r) : = $(e)(T)Q = e(0)Q + i^(s) [6{s) ® de{s) 

Jo 

+ (^Ric(X(s)) - it:(X(s)))ds]/i,t(M) (^(s)) ^ Q- (70) 

Let ^(r) be the right-continuous inverse of r — > 0(r); set D = {s : < 
s,A{s-) < A{s)}. If r > is fixed, then t(r) = A{(P{t)-) a.s. Now for g{T) an 
(^T-)-well measurable process such that r -E^[g'(r)^] is locally bounded then 
(see Theorem (6.6) in Ikeda and Watanabe [1]) 

E,[{r g{s)dB>'{s)}'] 

Jt{T) 

rA{u)As PT 

<E,[J2{ g{s)dB\s)f]=E,[ g{sfds]. 

udD ■^My--)^s Jo 

Then, for every T > there is a constant K = K{T) such that 

E,[\\m(r)\f]<K{l+ rE,m(s)\\']ds),yTe[0,T]. (71) 

Jo 

This proves that if ^ € S, then $(^) G S. Again, using Theorem 6.6 of [1], we 
have for ^i, ^2 G E, that 

EAmCi){r) - $(6)(r)f ] < k rEMi{s)-a{s)\f]ds,Te [0,T]. (72) 

Jo 
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Then, from standard arguments follows the following result. 

Theorem 3. The stochastic differential equation (60 — 63) has a solution 
K{t) € H which is unique. 

Let K{t) be the solution of (60 — 63); define 

M(t, w) := K{t, w) (g) e{T, w) = {K{t, w)%(t, w)), t > 0, (73) 

where we recall is the co-tetrad field dual to e. 

Theorem 4. M(r, w) is an i?" R^* ~ L{TM, i?"* ) -valued multiplicative 
horizontal functional of the generalized horizontal Brownian motion on Po{n) 
with reflecting boundary, i.e. 

(i)M(T, w) is (J^r) - adapted, (74) 

(M)Vr, s >0,M{t + s, w) = M{s, w)M{t, agw), almost surely (75) 

where the shift operator ag ■ M^(-fb(n)) ~^ ^(-fb(n)) is defined by {asw){T) = 
w{t + s) 

Proof: Part (i) is obvious by construction. We now fix s and consider 

K{t) := K(t + s, w),X{t) := X{t + s, w) and e(r) := e(r + s,w). Thus, K{t) 
satisfies eqts. (60 — 63) with respect to {X{T),e{T)). Still, if we apply the shift 
operator ag to K{t) we obtain that k{t) = K(t, agw) satisfies eqts. (60 — 63) 
with respect to (X(r, agii;), e(r, a^u;)) = (X(r),e(T)). If we set 

K'{t) := K{s, w) (g) 9{s, w) (g k{t), 

then, 

i^'(O) = K{s,w)®e{s,w)®{Ii^t^M){Xis,w))®e{s,w)^P + e{s,w)®Q) 
= K{s,w) ® {Iint{M){X{s,w)) ^ P + Q) = K{s,w), 

by (63) and the fact that ^ (g e = /, by definition. Hence K{t) and K'{t) satisfy 
both eqts. (60 — 63) and consequently, they coincide for every r. Then, 

K{t + s^w) = K{s, w) (8) 6{s, w) (g) K{t, agw), 

and thus multiplying by 9{t + s,w) = 0{T,asw) on the right yields the second 
part of this Theorem, c.q.d. 

An important property of M = (M{t,w)) is the following Lemma. 

Lemma 2. If X{0) e dM, then 

P®e{0)^ M{t) = 0, Vr > 0. (76) 
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Proof: It is enough to prove that P(g)6'(0)(8)iC(r) = 0,Vr > 0. IfX(O) G dM, 

then 

P0 0(O)0i^(O) = P®^(0)(7int(M)(^(0))e(0) ®P + e(0)®Q) 

= P(8)Q = PQ = 0. (77) 

Since K(r) = P(8>6'(0) (»K(r) satisfies eqts. (60-63) with respect to X(r), e(r)), 
then -fC(r) = by the uniqueness of solutions. 

Lemma 3. M{t,Taw) = AM{t,w)A-^,t > 0,yl G 0{n). 

Proof: Since X{t,Taw) = X{t,w) and e(r, T^tt;) = Ae{T,w), it follows 
from the uniqueness of the solutions of eqts. (60 — 63), that K{t,Taw) = 
AK{t,w). Thus, M{t,Taw) = AK{t,w) ^ 9{t,w)A-^ = AM{t,w)A-'^ . 

Proposition 1. For any r > such that X(r) G int(M), M(t) = M{t,w) 
satisfies 

dM(r) = M(r) ® [^Ric^ (r(r) - P^(r(r))]dr. (78) 

Proof: 

dM(T) = K{t) ® ^(r) [de{T) ® ^(r) + e(T) ^^(t) + de(T) • dOir)] 
+K{t) (8) ^(t)] ® [^Ric(r(r)) - P(r(T))]dr 

= ii:(r) ® 0(t) ® (i(e(r) 0(t)) + ii:(T) ® 0(T)[^Ric(r(r)) - R{r{T))]dT 

= K{t) ® e{T) ® [^Ric - i?](r(r))dT 

(since d(e(r)) (g) ^(r)) = (i(constant) = 0) 

= M(r) ® [^Ric^(r(T)) - P^(r(T))]dT, (79) 

where the last identity follows trivially from the definitions of Ric'' and E!' , since 
componentwise we have the equation for M = (M^) = {K^6^) 

dMi = {K%Rteie''^ + K%Rtleie''^)dT 
= {K%Rt5i + K%RfX)dT 
= {MlRi + M'^Rt,)dT, 
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where we finally note that the r.h.s. of (78) coincides with the coupling of M to 
the Weitzenbock term in which R has been substituted by E!' . 

Let Co(Po(n) ^ L{TM,R''*)) ~ T*M (g) r*Po(n) be the set of all bounded 
continuous functions F{r) on Po(n) taking values in L{TM, R"-* ), and such that 

P®9^F{r) = 0, if r = (x,e) G 9Po(n), ^ = e"\ (80) 

where again as before wc note that the (X)-product with F is meant in the image, 
E"* of the 2-form F(r) = ^F^adx"' A 6'"(a;) defined on TM A i?" by F. 
For oj e Co{Po(n) L{TM,R'^'')) and r > 0, set 

iHrLb){r) := Er[M{T, w) (g) uj{r{r, r, w))]. (81) 

Componentwise, we have that 

{Hr0j)ai3{r) := Er[M{T,wT^Cj^a{r{r,T,w))]. 

Theorem 5. {Fir} defines a one-parameter semigroup of operators on 
Co{Poin)^HTM,R^')). 

Proof: Prom Lemma (1) we obtain that for r G dPo{n), P®6®{HT0j){r) = 0. 
From the uniqueness of the solutions of the s.d.e. (60 — 63) we deduce the con- 
tinuity of the mapping r ^ € J-iW{Po{n)))j where J^{W{Po(n))) is the set of 
all probabilities on W{Po{n)) with the weak convergence topology; consequently, 
the functions {HT-Lj){r) are continuous in r. The equivariance of {Hr} follows 
from Lemmas 1 and 2. 

We are now in conditions of stating the solution of our boundary-initial- value 
problem (41&45), extending thus Theorem 1 to the case of smooth boundaries. 

Theorem 6. Let F{T,r) be a L(TM, i?"*)-valued smooth functions on 
[0, oo) X Po{n) such that for each r > 0, the mapping r F{t, r) is a function in 
Co{Po{n) L{TM,R''')). Then, with probability one, 

M(r) F(r, r(r)) - M(0) ® F(0, r(0)) = T M{s) ^ {L^F){s, r(s))(i5"(s) + 

Jo 

/ M{s) ® [— (s, r{s)) + HqF{s, r{s)) + (-Ric^(r(s)) - R\r{s)) ® F{s, r{s))]ds 

Jo OT 2 

+ M(5)®e(s)®Q®0(.)®[— (s,r(5))-^(s,r(.))r^^(X(5))e^(s) 

+ e{s) ® r„(X(s)) ® e{s) ® F{s, r(s))]#(s). (82) 
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Proof. The proof is an extension of a similar result for 1-forms, Theorem 6.5 
in [1], and is almost identical: One has to replace the expression ^Ric'' everywhere 
it appears by ^Ric'' — Ft . 

Note that the r.h.s. term of eq. (82) in which appears the coupling of 
F to the curvature terms , which consistent with our previous notations can be 
thought as an object of the form ^f^^^tix" coupling with ^Ric'' — i?'', coincides 
with the coupling of the Weitzenbock potential term ^Ric** — i? to an equivariant 
2-form F on Po{n) induced by F from the identity e F = F, or with our 
previous notation a 2-form Cl which satisfies = e<S)^^- Indeed, R^F = RF, and 
already Ric'' 0F = Ric^'-F, so that this term expresses the multiplicative operator 
functional for the zero-boundary case, while the last two terms in the identity 
(82) give precise account of the boundary conditions. Therefore, Theorem 6 can 
be regarded as a martingale problem solution [40] to the boundary-initial value 
problem (41&45); thus, we have proved that Ht-Co is the solution to this problem. 



4 The Navier-Stokes equations for the vorticity and 
the absolute boundary conditions 

We consider an oriented smooth connected manifold M with smooth bound- 
ary, dM, provided with a Riemannian metric, g, and with a time-dependant 1- 
form n(r, x) = Ur{x), r > 0, x G M. The Navier-Stokes equations for the velocity 
time- depend ant 1-form Ur{x) satisfying the incompressibility condition 

5ut-{x) = 0, 

is the non-linear diffusion equation Q 

du 

— (r, x) + VVl^Ur = lyVAiUr (83) 

where V denotes the projection operator into the co-closed component of the 
Helmholtz-Hodge decomposition of Ur (t > 0) which in view of the trivial identity 
(c.f. theorem 1.17 in [ 16 ]) 

Lu^Ur = V^gU^ + ^d{g{Ur,Ur))), 

^While in the case of a manifold without boundary, the viscosity term uAiUt {v is the 
kinematical viscosity) commutes with P, in the case with boundary we have to impose it from 
the beginning, i.e. we take i/PAiUr (c.f. page 144 in [17]) and in any case we are left with the 
expression in the r.h.s. of eq. (83) or eq. (84). 
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so that the non-hnear term is 

gUr = VL.fiUr, 

since by definition V vanishes on exact 1-forms, and thus the Navier-Stokes equa- 
tion takes the form 

du ^, . ^ , 

— = V{/^lUr + L_^9Ur) 

which from expression (10) we readily conclude that it can be rewritten as [4- 
6,12,13] 

-^{r,x) = VHi{2ug,-—Ur)ur{x),x eM,T>0. (84) 
We further have the no-slip boundary condition given by 

u\dM = 0. 

If we know define the vorticity 2-form of the fluid as 

Q, = flr{x) = dur{x), 

for any x € M, r > 0. Let us apply d to (84); in account that in writing down 
the Hodge decomposition of the non-linear term we obtain that 

and that similarly, further taking in account (11) we obtain that 

dVAiUr = A2nT, 

so that altogether we obtain the Navier-Stokes equations for the vorticity 

1 

— {T,x)=H2{2iyg,-—Ur)nr{x),xeM,0<T. (85) 

Since we have assumed that M has a smooth boundary dM, our interest now 
resides in boundary conditions. Then, 

n = ^napdx'' A dx^ 

where for simplicity we have ommited the variables (r, x), satisfies the absolute 
boundary conditions iff 

nan = 0,ya = l...n-l, and tt— O^fl = 0,1 < a < l3 <n. (86) 

ax" 
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These are the boundary conditions considered in Fluid Mechanics for the vorticity 
^ Now if we apply S to the definition of the vorticity, from eqt. (11) we obtain 
that 

5Qr = SdUr = — AlUr + dSur = —2Hi{g,0)ur, 

where in the last identity we have taken in account the incompressibility condition 
and eqt. (10) with Q = 0, and thus we obtain the Poisson-de Rham equation for 
the velocity 

Hi{g,0)ur{x) = -^6nr{x),T > 0,x G M. (87) 

We have already described in [5,6 & 13], the representations for this equation 
in the case of smooth boundary, by running either a gradient diffusion process 
by isometrically embedding M in an Euclidean space, or by running random 
Hessian and Ricci flows in the case of an arbitrary compact manifold, to solve for 
the vorticity and the velocity, respectively. Yet, in the above mentioned articles, 
the representation for the vorticity in the smooth boundary case was not given, 
so we shall now deal with this problem. 

We start by noting that from eqts. (1, 10&;85) and our constructions in 
Section I we can conclude that there is a Riemann-Cartan-Weyl connection on 
TM associated to the Navier-Stokes operator. This connection is determined by 
the Riemannian metric 2i^g, it is compatible with this metric, and has a trace- 
torsion determined by the time-dependant 1-form Q{t,x) = Qr{x) = ■^u(r, x) 
{x G M, r > 0) [4-6,13]. Thus, this "Navier-Stokes" connection, which we shall 
denote as V, has for Christoffel coefficients, the time-dependant functions on 
M defined by [| 

T%{t,x)=2v\^'^^{x) + ^-l-^|-i5^n(r,x)^ + i^5/37(^) «(^, } ,(88) 

so that we have a time-dependant horizontal canonical vector field, L^{T,r) = 
{L^ {t, r)), a = 1, . . . , n) associated to V by (4). It is trivial to check that expres- 
sion (88) indeed defines a connection compatible with the metric 2i'g and with 

^Indeed, from the no-slip boundary condition, it follows say, in the flat 3d case that the 
vector normal to dM given by Q.t ~ rot Ut (for every r) vanishes, which in invariant form is 
nothing else than Qan = 0, for every a = 1, . . . , n — 1. The condition (dfJT-)norm = 0, i-e.^ = 0, 
represents that the flow of vorticity through the normal to the boundary is zero, which is a 
natural condition for a non-permeable boundary, for every r. 

^The first term in (88) designate the Christoffel coefficients of the Levi-Civita connection V® 
defined in terms of the metric and its first order derivatives. 
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trace-torsion given by ^Ur- As long as eqt. (84) admits an unique solution in 
[0, T], then this connection is uniquely determined in [0,T] as well. 
We wish to solve the boundary-initial-value problem for 

neA\[0,T] xPo{n)),forT>0 

given by the Navier-Stokes equation for the vorticity as written on Po{n)- 

— (r, r) = H2{2ug, -_n,)17,(r), n{0, r) = Cj{r) (89) 

OT Iv 

where a) is the horizontal lift of a given 2-form a; defined on M (the vorticity at 
time 0), and [] 

e^(C^^a6 + C^)(T,r) =0, onaPo(n),Vr >0,a,6= l,...,n-l. (90) 
From (35) we note that we can rewrite the initial-value problem (89) as 

|-(r,r) = KiJ(r,r))(Lj)(r,r))(f).(r))+i.Ric^0l^(r,r) 

- 2z/^® fi(r,r),0 < T < T, (91) 

with 

fi(0)=tD, (92) 

and we already proved (c.f. equation (45)) that the boundary conditions admit 
the expression 

3Vt do!" 

p^e^n^ + Q^e^l-^-^^ r^^(r, x) ei{x) 

+ e(x) ®r„(T,x) 0(x) «)fi^](r,r) = 0,Vr = (x, e) € dPo(n), (93) 

where r„(r,x) = (T^p{T,x)), and (for < r < T) is in L{TM,R''*) and 
verifies that its partial horizontal lift e(x) ^rif) is the corresponding time- 
dependant 2-form on Po{n) iii the boundary-initial- value problem (91 — 93). The 

^These are the expressions on Po(n) of the given absolute boundary conditions which in p.d.e. 
theory are "Robin" -type conditions. 
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essential difference with (41&:45) is that the differential operator has a time- 
dependant drift term. 

Consider the Stratonovich stochastic differential equation for the diffusion 

2 

process r(r, r, t) = {X{x, r, t), e(e, r, t)) on Po{n) ~ >^ R"' defined by running 
backwards in time the s.d.e. (47); i.e. : 

dr^{r,T,t) = Ll{T-t,r{r,T,t))odB\t) 

+ 5"#(r - t), < t < r < T (94) 

(a = 1, . . . , n) with initial condition 

r(r, r, 0) = r, < t < r < T. (95) 

We note that it can be rewritten as 

dX'^{x, T, t) = e^(e, r, t) o dB''{t) + (5^#(r -t),0 <t <t <T (96) 

de^(e, T, t) = -rl^ir - t, X{x, r, t))el{e, r, t) o dX^ix, r, t) 
= -r^^(r - t, X{x, T, t))el{e, t, i)e(^ (e, r, t) o di?P(t) 

- r^^(r - t, X(x, r, t))el{e, r, t)#(T - t), < t < r < T, (97) 

(a, A; = 1, . . . , n) with initial conditions 

X{x, T, 0) = x, e"{e, r, 0) = e{x) G T^M, < t < r < T. (98) 

Here, </>(r) is a continuous non-decreasing process which increases only when 
X{x,T,t) G dM. |. Having assumed that g^p S C^{M), if we further assume 
that u{t,x) has all its components u{T,x)a G C^{M) (1 < a < n) and that 
further G C^(M), then (96 — 98) has a strong unique solution. 
For lD G Co(Po(n) ^ L(rM,i?"*)) and r G [0,r], set 

{HrCo){r) := Er[M{T,T,w) CdCo{r{r,T,T,w))]. (99) 

where M{T,t) = (M^(r, t, w)) is a i?" (8) i?" -valued process given by 

M(r, t, w) = K{t, t, w) (g) e{e, T,t),0 < t < t < T, (100) 
*Here we are assuming that B and (j> satisfy eqts. 55&56. 
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i.e. componentwise 

Ml{T,t,w) = KP^{T,t,w)e-p{T,t,w), 

where K{T,t,w) = {K^^{T,t,w)) is the unique solution of the s.d.e defined as 
fohows: Set K^{T,t,w) = K{T,t,w) P and K'^{T,t,w) = K{T,t,w) Q so 
that K{T,t,w) = K^{T,t,w) + K'^{T,t,w), where for any r G [0,T] such that 
X{x,T,t) G int(M), t £ [0, r] we have 

(i) dK^ (r, t) : = dK{T,t) ® P = K{t, t) ® [9{e, r, t) O de{e, r, t) 

+ {i^Ric(X(x, r, t)) - 2vR{X{x, r, t))}c/t] ® P, (101) 

and for any r € [0, T], 

(ii) dK^{T,t) : = dK^{T,t)®Q 

= K{T,t)^[9{e,T,t)(8)de{e,T,t) 

+ {i/Ric(X(x, r, t)) - 2uR{X{x, r, t))}dt]/int(M) (^(a^, r, t)) (g) Q, (102) 

and with probabihty 1, r ^ X^(r, t) is right-continuous with left-hand limits: 
furthermore 

K^(t, t) = 0, if X{x, r, t) G {x" = 0}, (103) 

and the initial value is 

Ki(r,0) =/i,t ^M){X{x,T,0))e{e,T,0) ^ P,K\t,0) = e{e,T,0) (^Q. (104) 

Theorem 7. Let a) = (I'(r) a be function in Co{Po(n) L{TM,R^*)). 
Then, (ifT-tD)(r) = Er[M{T,T,w) u;(r(r, r, ^i;))] is the unique solution of (91), 
with initial value 0(0, — ) = u;(— ), and further satisfying the boundary condition 
given by (93). 

Proof: It follows easily from observing that this is the time-dependant 
version of Theorem 5. 

Now, let $(r, x) = $r(a;) G A2([0, T]y.T*M) such that its horizontal (partial) 
lift ^t{x) e{x) = {HT-io){r), for all r = (x,e(x)), O < r < T. Then, from 
Theorem 7 follows that $ is the unique solution of the boundary-initial-value 
problem 

— {T,x)=H2{2ug,^Ur)^rix), ^oix) = uj{x) ,0 < T < 0, X £ M , (105) 

^We shall omit in the following the variable w, as customary, in the expression for K and its 
specialization to the flat case. 
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where u) is the 2-form on M such that its horizontal lift is a), i.e. u>{x) (8) e(x) 
e{x) = u}{r), and satisfying further the absolute boundary conditions: 

^an = 0, -^^ap = 0, 1 < a < (3 < n - 1. (106) 

Thus, we have obtained the representation for the Navier-Stokes equations for 
the vorticity with absolute boundary conditions. 

Now we can solve for the vorticity equation in the flat space -R" . Now since 
in this case g = I, the identity, and thus in (88) the expression for the Christoffel 
symbols reduce to the form 

r^7(^'^) = {^^P «(^'^)7 + ^^/37 <r,xr^ , (107) 

so we now consider the s.d.e. 

dX'^ix, r, t) = e^(e, r, t) o dB''{t) + (5^#(r -t),0<t<T, (108) 

det{e,T,t) = =-T'^^{T-t,X{x,T,t))el{e,T,t)e^p{r,T,t)odBP{t) 

- r^^ir -t,X{x,T,t))el{e,T,t)dct>{T -t), (109) 

{a,k = 1, . . . ,n) with initial conditions 

X{x, T, 0) = X, e"(e, r, 0) = e{x) G T^M, 0<t<T. (110) 

Here, 4>{t) is a continuous non-decreasing process which increases only when 
X(r, t, x) € {x" = 0}. Furthermore, since Ric and R are identically equal to in 
the expressions (101&102), so that K(T,t) = K^{T,t) + K'^{T,t) is the solution 
of the problem described as follows: for any r G [0,r],x G dM we have 

dK^{T, t) := dK{T, t) P = K{t, t) (g) 6'(t, t) de{T, t) (g) P, (111) 
and for any r G [0,T], 

dK^{T, t) := K{t, t) (g) ^(r, t) ® de{T, t)/i„t(M) (^(r, t)) Q, (112) 
satisfying 

K^{T,t) = 0,ii X{T,t) & {x"" = 0}, (113) 
and the initial value is 

K\t, 0) = /i,t(M) (^(r, 0))e(r, 0) g) P, i^2(r, 0) = e(r, 0) Q. (114) 



28 



Thus, if we finally consider M{T,t,w) = K{T,t,w) (8) 9{T,t,w) (where = e ^) 
for u e Co{Poin) L^BJ", i?"*)) and < r < T, set 

{HrUj){r) := Er[M{T,T,w) 0u>{r{r,T,T,w))]. (115) 

Theorem 8. The 2-form G A^([0,r] x i?"*) such that 

^t{x) (8) e(x) = {Hroj){r), r = {x, e{x)), (116) 

is the unique solution of the boundary-initial-value problem given by 

1 

— = H2{2vl, ——Ut)^t = i^div grad fir — Lu^Vlr (117) 
with initial value 

e(0) (8) $(0, -) = a;(-) (118) 
and further satisfying the absolute boundary conditions . 

5 Kinematic Dynamo of Magnetohydrodynamics on 
Smooth Manifolds with Smooth boundary 

We can extend these constructions to solve the passive transport equations 
of magnetohydrodynamics. All we have to do is to construct the multiplicative 
horizontal operator functional for the differential operator Hn-i{2i/"^g, i^^Ur), 
with z^™" the magnetic diffusivity, Ur given by the solution of cither the Euler 
equations, or still the Navier-Stokes equations, with absolute boundary condi- 
tions. Indeed, doing this (for which we have to take in account the coupling of 
the Ricmannian curvature to the (n — l)-magnetic form in the Weitzenbock term 
for the definition of K), we can then solve the kinematical dynamo problem for 
a magnetic field B{t,x) defined by duality as 

is,// = 0^,0 < T < r (119) 

where i denotes the interior product derivation, ^ = det(g)2c?x^ A ... A dx" is 
the Riemannian volume density, and Q{t,x) G A"~^([0,r] x T*M) satisfies the 
transport equations (kinematic dynamo problem) 

dVl -1 

— = Hn-l{2v'^g, -—Ur)nr,no = iBoH,0 < T < T, (120) 
OT Zu^ 

for given 

B(0,-) = B(-), (121) 
satisfying the absolute boundary conditions for Q. 
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6 Conclusions 



In this article, we have given random representations for the Navier- Stokes equa- 
tions for the vorticity, and the kinematic dynamo equation of magnetohydrody- 
mamics, for incompressible fluids on smooth boundary manifolds, and in partic- 
ular, in the case of flat cuclidcan space. No such general representations were 
known but for the case of empty boundary smooth compact manifolds [5,12,13], 
and as an implementation of the case of manifolds isometrically immersed in 
Euclidean space, for Euclidean space itself. We would like to remark that this 
program stemmed as a covariant extension of the random vortex method in 2D 
(without boundary) of Computational Fluid Dynamics ( see A. Chorin [32], and 
references therein). In the latter method,due to the fact that in fiat 2D the vor- 
ticity can be identified with a scalar field, the vorticity equation can be integrated 
by using the Ito formula for scalar fields. The extension of this formula, has been 
the backbone for the obtention of our representations in both the empty and 
non-empty boundary cases. We should also stress that product formulas (such 
as those arising from transition densities as is the case for diffusion equations) 
have been algorithmically implemented in Fluid Dynamics for a long period [33]. 
Yet, more specifically related to the present approach, numerical methods for the 
random integration of nonlinear partial differential equations have been devel- 
oped (c.f. [34]); furthermore these methods have been implemented for p.d.e.'s 
for scalar fields satisfying reflecting boundary conditions [35-38]. Thus we might 
expect that the former approach properly extended to differential forms and im- 
plemented for the representations achieved in this article, will eventually lead to 
interesting numerics. 
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